RlC'd PCT/PTO 0 8 AUG 2005 



DE VICE AND SOFTWARE PACKAGE FOR EXTRACTING A GEOLOGICAL 
HORIZON AND RELATED PROPERTIES 



The invention relates to a process to extract a 
geological horizon and related properties from their image 
deduced from seismic data. 

The invention also relates to a device to extract a 
geological horizon and associated properties from their 
image deduced from seismic data. 

The invention finally relates to computer software to 
extract a geological horizon and related properties from 
their image deduced from seismic data. 

There are. known processes for extracting the horizon 
and their seismic attributes. In these known processes, 
there ■ is constructed a three-dimensional seismic 
attribution matrix of the same size and same sampling as 
the initial three-dimensional seismic matrix. 

To attribute to each point on the horizon a value of 
seismic attribution, one can either interpolate the 
surrounding modes of the three-dimensional seismic 
attribution matrix or select the stored value in the mode 
which is closest to the three-dimensional seismic 
attribution matrix. 

In the case of interpolation, the corresponding 
preprocessing requires a large processing time and a very 
great quantity of available memory; this is why it is 
generally necessary to carry out this pretreatment on a 
powerful computer. 

In the case of the choice of attribution of the 
nearest node, this pretreatment is not necessary, and the 
determination of the attribution can only be carried out at 
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points adjacent the points of the network of the original 
three-dimensional seismic matrix. 

The principal drawback ' of the prior art is to 
introduce . the defects or artifacts which induce critical 
5 errors on the horizon and the seismic attributes and lead 
to .poor geological interpretation of the seismic 
measurements . 

The first object of the invention is to improve the 
known technique by providing a direct and more exact 

10 computation, by minimizing the risk of possible error on 
the attributes at all points on the horizon. 

A second object of the invention is to suppress the 
vertical errors that are adapted to lead to a poor 
geological interpretation . 

15 The invention has for its object a process to extract 

a geological horizon and associated properties from a 
seismic representation, in which there is constructed a 
continuous function S ij/k (t) by interpolation or 
approximation of the discrete seismic traces of a 

20 multidimensional seismic matrix, said function being 
designated as a "continuous local seismic trace", 
comprising the following steps: 

a) . using as optimum offset of two adjacent continuous 
local seismic traces, the value of the offset rendering a 

25 maximum their correlation function; 

b) . taking as conditional neighborhood of a continuous 
"central" local seismic trace S ljfk (t) the sub-neighborhood 
consisting in adjacent traces S pq , k (t) corresponding to 
optimum offsets associated with correlations Rij fP q f k(h) 

30 greater than a predetermined threshold comprised between 0 
and 1; 
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c) . constructing a two-dimensional extraction matrix 
adapted to be filled with extracted values; 

d) . choosing a seed point P(i,j,t) and determining the 
point P(i,j,k) that is vertically closest; 

5 e) estimating the related properties of the 

conditional neighborhood and filling the two-dimensional 
extraction matrix with offset properties by translation of 
the current variable (t) of the optimum offset value 
(hij, pq#k ) corresponding to the point P(i,j,k) vertically 
10 nearest. 

According to other characteristics of the invention: 

- the content of the two-dimensional matrix of step c) 
is successively replaced in the course of successive 
repetitions of extraction, 

15 - at each repetition, there is utilized as new seed 

points of step d) all the points corresponding to locations 
of the two-dimensional matrix of step c) filled with the 
preceding repetition, 

- the property of extracted subsurface is the seismic 
20 amplitude reflected and detected by geophones, 

there is visualized on a visualization screen 
horizon surfaces painted with colors corresponding to a 
code of the extracted properties. 

The invention also relates to a device for practicing- 

25 the process according to the invention, comprising means to 
use as optimum offset of two adjacent continuous local 
seismic traces, the value of the offset rendering a maximum 
their correlation function, means to take as conditional 
neighborhood of a central continuous local seismic trace 

30. Sij, k (t) the sub-neighborhood consisting in adjacent traces 
Spqk(t) corresponding to optimum offsets associated with 
correlations Ri jfP q f k(h) greater than a predetermined 
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threshold comprised between 0 and 1, means to construct a 
two-dimensional extraction matrix adapted to be filled with 
extracted values, means to select a seed point P(i,j,t) and 
to determine the point P(i,j,k) which is vertically 
5 nearest, and means to estimate the related properties . of 
the conditional neighborhood and to fill the two- 
dimensional extraction matrix with offset properties by 
translation of the current variable (t) of the value of the 
optimum offset corresponding to the point P(i,j,k) which is 
10 vertically nearest. 

According to other characteristics of the invention: 

- the device comprises means for memorizing and means 
for visualizing seismic parameters determined with the help 
of the process according to the invention. 

15 The invention finally has for its object a computer 

software, comprising elements of program code to carry out 
the steps of a process according to the invention, when 
said program is executed by a computer. 

The invention will be better understood from the 

20 description which follows, given by way of non-limiting 
example, with reference to the accompanying drawings, in 
which : 

- Figure 1 is a schematic representation of a three- 
dimensional seismic matrix and a continuous local seismic 

25 trace. 

- Figure 2 is a schematic representation of an example 
of local neighborhood N ijfk consisting of an assembly of 
continuous local seismic traces S pq , k (t)' adjacent the central 
reference trace S ijfk (t) itself located on the vertical of 

30 the geophone G ±j and interpolating the seismic amplitudes 
about t=k . 
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Figure 3 shows schematically a vertical cross- 
section of conditional local neighborhood N ij/k (r) of a 
continuous reference local seismic trace S ij/k (t). 

Figure 4 shows schematically a vertical cross- 
section of a three-dimensional seismic matrix with optimum 
offset hij, pq , k and a conditional neighborhood N ij/k (r). 

Figure 5 shows schematically a vertical cross- 
section of a three-dimensional seismic matrix analogous to 
Figure 4 with the presence of a fault or discontinuity. 

- Figure 6 shows schematically a functional organogram 
of a process according to the invention. 

- Figure 7 shows schematically an enlarged partial 
view of a three-dimensional seismic matrix having a 
discontinuity and illustrating the practice of the process 

15 according to the invention. 

With reference to Figures 1 to 5, identical or 
functionally eguivalent elements are designated or marked 
in an identical manner. 

In Figure 1, a three-dimensional seismic matrix is 
obtained by picking up the measurements registered by the 
geophones G ±j disposed on a network x,y at coordinate points 
i, j. The pickup of these samples measurements time-wise is 
represented along a descending axis t representative of the 
depth of a vertical descending from the surface of the 
ground or of the sea. The measurements are characterized 
by their amplitude, for example an amplitude picked up by 
the geophone G ±j with time or at the sampling depth t k . The 
discrete measurement carried out by the geophone Gj.j in time 
or at depth t k is called the seismic amplitude S ij/k . 

The assembly of seismic amplitudes corresponding to a 
geophone G ±j of coordinates i, j is a one-dimensional matrix 
(Sijl, Sij2, S ijk , S ijN ) called a discrete seismic 
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trace, because this one-dimensional matrix corresponds to 
the trace according to the point of horizontal coordinates 
i, j of the three-dimensional seismic matrix obtained by 
seismic measurements. 
5 The vertical axis t oriented along a descending 

vertical designates usually the time, but can also be 
considered to represent a dept from the surface. 

The invention relates not only to the application of a 
third coordinate t representative of time, but also to a 

10 third coordinate t representative of depth. 

From the discrete seismic trace located on the 
vertical of a geophone G ±j , there is defined, by 
interpolation or approximation, discrete values about 
t _t k =k, a continuous function S ij/k (t) which is designated as 

15 the "continuous local seismic trace". The approximation or 
interpolation methods for discrete values to give rise to a 
continuous function are numerous, and comprise particularly 
polynomial interpolations or approximations, as well as 
polynomial trigonometric interpolations or approximations. 

20 Any other interpolation or approximation variant 

giving a continuous function can also be applied to the 
present invention to provide a "continuous local seismic 
trace" . 

In Figure 2, several continuous local seismic traces 
25 define a neighborhood of a central reference continuous 
local seismic trace S ij/k (t). The neighborhood of a 

continuous local seismic trace S ij#k (t) is defined as the 
assembly of the continuous local seismic traces whose 
horizontal spatial indices p, q are adjacent horizontal 
30 spatial indices i, j of the reference continuous local 
seismic trace. 
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By way of example, the horizontal spatial coordinates 
p, q corresponding to geophone G pq are adjacent the 
horizontal spatial coordinates i, j corresponding to the 
geophone G±j if the absolute values of the differences i-p 
and j-q are less than given whole numbers, for example 2. 

In this case, as shown in Figure 2, the continuous 
local seismic trace S ij/k (t) is associated with eight 
neighboring continuous local seismic traces surrounding the 
central continuous local trace Sij, k (t). 

In the case of continuous local seismic traces 
produced by seismic measurements, the profile of the 
geological horizons introduces vertical offsets between the 
adjacent continuous local seismic traces.. For determining 
the relationships between two adjacent local continuous 
seismic, traces S lj#k (t) and S pq , k (t) centered on. the same 
sampling vertical coordinate t=t k =k and corresponding to 
different spatial coordinates i, .j and p, q, there is 
calculated the correlation function R ij/Pq , k (h) of two 
continuous local seismic traces. 

The correlation function of two adjacent continuous 
local seismic traces is obtained by the following equation 



in which the numerator corresponds to the covariance 
function of S ij/k (t) and S pq , k (t) obtained by the following 

expression <W<h) - JS M (t). S^t + h).*. 

flc-A 



in this integral defining C ij(Pq , k (h) , the parameter A defines 
a "vertical investigation window" about t=t k =k. 

For example, if the continuous local seismic traces 
S i3(k (t) and S pq , k (t) are trigonometric polynomials of the 
following form interpolating the seismic data 



S^t) - £a«* .cos(s»t) + b. OJt .siii(sot) 
%k(t) - jTa^cos(s©t) + b^sin(sa>t) 



then one can show that the covariance function d^pq.ieth) is 
itself a trigonometric polynomial of the following form in 
which the coefficients A«'«' k and B«'«' k depend on the 
coefficients a s ^' k , a s ^' k , b s ^< k and b s ^ k : 



A conventional and known mathematical result is that the 
correlation function 



•o n . C a,r»* (h) translates a similarity of the 

correlative functions when this correlation function 
approaches 1 . 

The study of the correlation functions of the 
continuous local seismic traces S ij(k (t) and S pq , k (t) permit 
defining- as optimum offset the value h ij<pq , k of h 
corresponding to the maximum of the correlation function 
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Rij, pq ,k(h), which is to say to the maximum of the correlation 
function nearest 1. Contrary to the prior art based on a 
discrete formulation of the correlation function Rij, pq , n (h) , 
the use of a continuous formulation of the seismic traces 
5 and hence of the correlation function permits obtaining an 
optimum offset which is not necessarily a whole number 
multiple of the sampling interval along the vertical axis 
corresponding to the variable t. 

The definition of the optimum offset h ij(Pq , k of two 
10 adjacent continuous local seismic traces permits obtaining 
a first approximation of the horizons passing through this 
reference continuous local seismic trace S ij(k (t), as 
represented in Figures 3 to 5. Contrary to the prior art, 
the fact of using optimum offsets which are not necessarily 
1.5 a whole number multiple of the sampling interval along the 
vertical axis, permits avoiding errors known as "aliasing". 

In Figures 3 to 5, there is defined a conditional 
neighborhood N ijfk (r) of the reference local continuous 
seismic trace S ij(k (t) which is a sub-neighborhood of the 
20 initial neighborhood of the continuous local seismic traces 
used for the computation of correlation and definition of 
the optimum offsets. The conditional neighborhood N ij<k (r) 
is selected such that, for any local continuous seismic 
trace S pq , k (t) belonging to- N ij(k (r), the corresponding 
25 optimum offset h ij(Pq , k is such that the correlation 
Ri j . P ,. k (hi j . pq . k ) between S ij;k (t) and S pq , k .(t) is greater than a 
predetermined threshold r comprised between 0 and 1. 

As can be seen in Figure 5, the fixing of the 
threshold r also permits encompassing in a conditional 
30 neighborhood a discontinuity or geological fault, which 
constitutes an important advantage relative to the prior 
art . 
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The above arrangements thus provide a continuous 
modeling permitting the practice of the invention. 

With reference to Figure 6, a process according to the 
invention comprises a step of digital modeling 100 to 
define continuous local seismic traces, optimum offsets and 
conditional neighborhoods, as described with reference to 
Figures 1 to 5. 

In step 101, there is defined a two-dimensional matrix 
whose indices of lines and columns correspond to the 
coordinates of the seismic * measurement geophones. 

In step 102, there is selected a seed' point P(i,j,t) 
or first extraction point. 

In step 103, there is determined the node (i,j,k).of 
the three-dimensional seismic matrix nearest the seed point 
P(i,j,t). 

This node has the spatial coordinates i,j and time or 
depth coordinates the coordinates t k =k nearest the 
coordinate t of the seed point P(i,j,t) selected in step 
102. 

In step 104, there is selected first of all the local 
continuous local seismic traces S pq , k (t) belonging to the 
conditional neighborhood N ij/k (r) defined in step 100 and 
such that the indices (p,q) correspond to an empty location 
on the two-dimensional matrix defined in step 101. 

Then, for each local continuous seismic trace S pq , k (t) 
thus selected, there is attributed to the point P(p,q,t) 
the value P(p,q,t + h ij/Pq , k ), in which h ij/Pq , k is the optimum 
offset of the point P(p,q,t + hij,pq,k) relative to the 
adjacent point P(i,j,k) determined in step 103. This point 
P(p,q,t + hij, pq , k ) is considered as a new point located on 
the horizon and is stored at the location (p,q) on the two- 
dimensional matrix defined in step 101. 
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The process which has been described is adapted to 
determine the assembly of the points belonging to the 
horizon passing through a point P(i,j,t), but it also 
applicable to attribute to these points, to the extent of 
5 their determination, an assembly of properties 
characterizing the physical nature of the terrains adjacent 
the horizon thus determined. 

To this end, one can attach to each point 
P(p,q,t + hij #P q, k ) of the two-dimensional matrix defined in 
10 step 101, an assembly of properties called ""seismic 
attributes' 7 . 

These properties are generally represented by 
different colors on screen displays, so as to permit a 
rapid visualization without risk of error of the properties 

15 associated with a geological horizon. 

In practice, there are - carried out steps similar to 
the steps described with reference to steps 100 to 104 and 
there are constructed the attributes of the point 
P(p,q,t + hij,p q , k ) with the help of the equation of the local 

20 continuous seismic trace S pq , k (t); for example, there is thus 
constructed an envelope attribute, a phase attribute or a 
frequency attribute, in a manner known per se, for example 
according to information in the publication of TANER M . T . , 
KOEHLER F. , SHERIFF R.E., (1979), Complex seismic trace 

25 analysis, GEOPHYSICS, volume 44, No. 6, pages 1041 to 1063. 

This method thus permits displaying the seismic 
attribute on the extracted horizon by painting ■ it with 
corresponding colors . 

The process which has been described thus permits 

30 defining the extraction value's of a predetermined horizon 
corresponding to real data, but only in a conditional 
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neighborhood of the point P(i,j,k) nearest the selected 
seed point P(i,j,t). 

To propagate the process of extraction to the assembly 
of two-dimensional matrix defined in step 101, there is 
first of all constructed an assembly of new seed points 
constituted by new points P(p,q,t + h ij(Pq/ ,.) located on the 
horizon which have been determined and are stored in 
corresponding locations (p,q) of the two-dimensional matrix 
defined in step 101. Each of these points is memorized in 
step 106, to be re-injected in step 102, and there serving 
as a new seed point. 

This propagates the successive processing to the 
assembly of the two-dimensional matrix corresponding to the 
geophones . 

When no new seed point can be found, the process is 
oriented in step 106 toward a step 107 of visualization on 
the screen of the . horizons or attributes or properties 
associated with the geological horizon, to as to permit a 
visualization in color or in gray of the corresponding 
properties of the corresponding extracted geological 
horizon corresponding in real time t. 

Thus, at the end of the algorithm defined in Figure 6, 
the two-dimensional matrix of step 101, contains a 
plurality of points corresponding to the same horizon as 
the selected initial seed point P(i,j,t). 

The assembly of the points can be reordered with its 
neighbors, according to Figure 7, to define a network 
entirely covering the horizon containing the seed point 
P(i,j,t), this network thereby permitting representing the 
extracted horizon as a locally continuous surface on the 
visualization screen of a work station. 
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The process according to the invention is preferably 
practiced with a device comprising suitable means for 
practicing the steps described with reference to Figure 6. 

In particular, a device according to the invention 
comprises memory means needed for the successive 
computations and the visualization means necessary for step 
107. 

Preferably, a device according to the invention is a 
programmable device controlled by computer software 
comprising program code elements to execute the steps of 
the process described with reference to Figure 6. 

The invention described with reference to several 
particular objects is in no way thereby limited, but on the 
contrary covers all modification of form and any variation 
of embodiment within the spirit and scope of the invention. 
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